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STAP-8; SPIKE TRAIN ANALYSIS PROGRAM 


DECUS Program Library Write-up 


DEC US NO. 8-342 


1. ABSTRACT 

STAP-8 is a subroutine package for basic statistical analysis of stochastic point processes, 
written in PDP-8's PAL symbolic assembler language and available as a machine language 
perforated tape in binary format. Main effort was not made on sophisticated statistical 
techniques, but rather on a easy to control, variable program library, including CRT- 
display, paper tape output and teletype listings. 

2. REQUIREMENTS 

Storage - The package occupies all of the basic 4K core memory's locations, except the 
last page (BINLOADEfy^MONITOR). During load procedure, the package's analysis pro¬ 
grams are saved on a full DF32 random access disk file track, service and I/O routines only 
rest in core permanently. Saved core space will then be loaded with DEC'S basic floating 
point package and again saved on the lower half of a second DF32 track. (For detailed 
storage information, see MEMORY CHART appendix.) 

Equipment - Basic PDP-8, 4K core memory;EAE extended arithmetic element;DF32 random 
access disk file (3 tracks); A>08 laboratory peripheral; CRT display unit. 

3. USAGE 

3.1 Loading - The package may be placed in memory either by the use of the Binary Loader 
or the system monitor's equivalent .LOAD. 

NOTE: - Throughout this STAP-8 manual, machine language memory addresses and instruction 

codes are written in a four-character octal format '. 

The user with little or no experience may therefore make use of the octal-binary conversion 
table, listed below. 

0 000 4 100 

1 001 5 101 

2 010 6 110 

3 011 7 111 

Now, to load STAP-8 proceed step by step as fol lows: 

.1 Set PDP-8 front panel switch register (SR) to 7600 , then depress LOAD/ADD and START. 

If the teletype outputs a ". ", continue as in step .2'. Else, consult the ERROR DIAGNOSTIC 
part of this manual 1 . 











.2 Type 


.LOAD 

* IN—T: 


* 


*OPT-2 


ST=34j0(0 


NOTE: - Each line should be terminated with a RETURN command; wait until monitor types 
the next option request code l 

.3 LOAD responds then, typing an upward arrow. Turn low-speed reader control (LSR) to 
position FREE and put STAP-8 paper tape in the reader. 

.4 Type CTRL/P and turn LSR control to START. The loader will now read the tape and 
halts after checking the final checksum. If no checksum halt occurs, i.e., if the reader 
continues until the tape end has passed through, consult ERROR DIAGNOSTIC parti 

.5 Turn LSR to FREE, reload STAP-8 paper tape, depress CONT, type CTRL/P and again 
turn LSR to START. 

.6 .LOAD performs second-pass loading and again shows halt after checksum diagnostic. 

.7 Turn LSR to FREE, depress CONT and type CTRL/P. 

.8 Program control now turns over to STAP-8 disk save subroutine, which introduces itself 
by the following character string output: 

SET SR FOR LOADING MODE / THEN CONTINUE 

.9 Set switch register to 7600 and depress CONT. STAP-8 analysis programs will be saved 
now on a full DF32 disk track, and, after transfer complete control returns to the system's 
monitor, which again responds with a ". ". 

.10 Proceed now as in step .2, but replace the restarting address by: 




ST=3447 


.11 Procedures for two-pass loading the DEC 8-5-S/A or equivalent DEC 8-25-F/A floating 


point package are the same as described in step 3-7. 

.12 After typing the last CTRL/P, control turns over permanently to STAP-8; string output 
announces: 


PROGRAM LOADED / SELECT SUBROUTINE: 
TYPE AND CODE: 


and lists the mnemonics of the package's 11 analysis programs available! 
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3.2 BINLOADER mode 


If the system's monitor fails to work (see ERROR DIAGNOSTIC), STAP-8 may be loaded by 
the use of the RIA^/BJN loader; proceed as in the steps listed below. 

.1 * Toggle in the RIM loader, using the front panel switch register and the instruction 
codes listed on the left side of the panel l 

.2 Set SR to 7756, depress LOAD ADD and START. Put the BINLOADER paper tape in the 
reader and turn LSR control to START. 

.3 After the reader halts at the end of tape (no checksum diagnostic occurs 1 .) depress STOP, 
set SR to 7777, depress LOAD ADD and START. 

.4 Put the STAP-8 paper tape in the reader and turn LSR control to START. The BINLOADER 
places the program in the memory during one-pass only; a checksum halt should occur before 
end of tape has passed through! 

.5 Set SR to 34,0,0, depress LOAD ADD and START; the teletype then should print the message 
mentioned in step 3.1.8 

.6 Set SR to 7777 and depress CONT and then put the DEC 8-5-S/A or DEC 8-25-F/A 
floating point package in the reader, turn LSR control to START. 

.7 After one-pass loading and checksum halt, set SR to 3447, depress LOAD ADD and 
START. The teletype output will be the same as mentioned in step 3.1.12 before! 

3.3 - Operating Procedures 

3.3.1 Program Flow 

There is no HLT (halt) instruction in any part of the program, so the user will be discharged 
of front panel control operations, except for AX08 interface connections. The only way to 
control the program flow is realized by the mean of teletype 'shorthand' intercommunication . 
This will avoid hazardous destroying of program instruction by the user, setting incorrect 
switch register or control keys. If the RUN indicator on the panel is turned off, the program 
control has executed a HLT instruction, and the STAP-8 program flow is broken . 

NOTE: - NEVER TRY TO RESTART THE PROGRAM, BUT RELOAD IT AS DESCRIBED IN 3.11! 

3.3.2 Program States 

The program flow may be divided in three states with respect to the teletype intercommunication. 

OLS - (Off-line state) Teletype keyboard information is not sensed during program evaluation 
and therefore will be ineffective on program flow. 


* 


For details, see description in PDP-8's software manual! 











PRS - (Program request state) Program control turns over to the teletype input subroutine 
and expects the user to type a message. Control will not return until message terminating 
RETURN has been typed'. 

URS - (User request state) During program evaluation, the teletype input flag periodi¬ 
cally will be sensed for auser request, which will cause a jump to the teletype input 
routine (see PRS). 

3.3.3 Teletype Intercommunication 

At the end of loading procedure, i.e., after instruction code mnemonics listing com¬ 
pleted, control turns to PRS (see 3.3.2) and expects the user to type: 

- a leading "*" , 

- one of the three-character instruction codes, 

- and the terminating RETURN keyl 

This calling sequence will load the corresponding program from its save area on DF32 disk 
into core memory locations 3400 - 3577, and program evaluation continues immediately 
at 3400. The program will introduce with an identifier label and continue with comment 
or request printout, or start data processing. 

This is the general rule how to monitor program flow control; further details are mentioned 
in section PROGRAM DESCRIPTION. 

Also, error messages may occur, discussed in the section ERROR DIAGNOSTIC. 

4. PROGRAM DESCRIPTION 

STAP-8 is a collection of subroutines; a part of it, composed of input-output, CRT 
display, format-conversion and other service programs may not be controlled or used 
through teletype intercom command mode, and therefore are not discussed below. 

The 11 subroutines, referred to as 1 the program 1 which compose the 'core 1 of STAP-8 
package and compute some of the elementary statistical measures of interspike time interval 
series, are described and discussed in detail. 

Although the programs were written for neuronal 'spike' potential analysis, and therefore 
make use of the versatile AX08 laboratory peripheral, a STAP-8 version will be available, 
where raw data input consists of any digital data files up to 1,024 12-bit words each. A 
data word will then be interpreted as a positive integer in the range of 0 to 4095, repre¬ 
senting a time interval between two successive events of a series, considered to be the 
realization of a stochastic point process'. 


4.1 OPT - OPTions 

This program includes output and service subroutines for various analysis extensions. 






OPT introduces with label and colon printout, and expects the user to select and type 
(PRS) one of the legal option codes listed below, and a terminating RETURN 1 . 

RC - AXJ08 RC-clock time base calibration. 

ST - AX08 Schmitt-Trigger input level control. 

DN - Delete old N and accept new interval file range. 

P - Punch binary coded paper tape output. 

L - List statistical measures of central tendency. 

Non-legal codes will again cause the label and colon printout, i.e., program restart! 

If the code has a leading the command decoder interprets this as a user's program 
request; OPT will then be left and a new program loaded! 

4.1.1 RC 

Code terminating RETURN initiates the program to set AX08 digital output S0 from 0 to -3 V 
at the AX08 RC-clock pulse rate. Link S0 output to an oscilloscope, set the vertical 
amplifier to 1 V/cm and the time base generator to the range desired. The time base is 
equal to the duration of one positive or negative going square-wave pulse. Calibrate the 
square-wave by the use of RC-clock control potentiometers on the AXJ08 front panel! 

NOTE: - TIME BASE CALIBRATION REQUIRES CAREFUL ATTENTION. THE RESULTS 
OF STATISTICAL COMPUTATIONS ARE RELATIVE AND BASED UPON THE TIME UNIT 
SELECTED! 

4.1.2 ST 

Code terminating RETURN initiates the program to type the input channel request "S:", 
which the user has to respond, typing the Schmitt-Trigger input number (1, 2 or 3) he wants 
to test, and terminating RETURN. 

The program then senses the state of the Schmitt-Trigger selected, and, if set to 1, pro¬ 
duces a -3 V pulse (0.5 msec duration) at S0. S0 may be linked to an oscilloscope and, 
if there is a dual beam option, the original spike train may be compared with that, recorded 
by the Schmitt-Trigger. 

Schmitt-Trigger level will be set from 0 to -2.5 V by Schmitt-Trigger potentiometers on 
the AXJ08 front panel! 

NOTE: - BOTH RC AND ST OPTIONS MAY BE INTERRUPTED, STRIKING ANY KEY ON 
THE TELETYPE KEYBOARD. THIS WILL CAUSE A NEW LABEL AND COLON PRINTOUT! 

4.1.3 DN 

Code terminating RETURN initiates the program to type the new file range request "N:", 
which the user will have to respond, typing the upper boundary of the next spike train to 
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be processed and terminating RETURN 1 . 

Core location ,0,022 will then be loaded with the new N's two-complement value; the 
value of the old N is lost 1 . 


NOTE: - N SHOULD BE LESS THAN 1024'. 
4.1.4 P 


Code terminating RETURN initiates program request typeout "I:". The user is expected 
to design the data file in core memory to be punched on paper tape in binary format. 


I:R Output source is the 128x3 word register. It contains the integer (floating point) results 
oFstatistical computations. The memory chart below and program descriptions of the 
corresponding statistical measures will help the user to plan additional off-line analysis 
programs. 


REGISTER MEMORY CHART 


2600 2601 2602 2603 2604 2605 3375 3376 3377 



Serial correlation coefficients (SCG) will be computed and stored in floating point format: 


EXP 11-bit exponent plus sign bit 

HORD High-order 11-bit mantissa plus sign bit 

LORD Low-order 12-bit mantissa 


INTEGER 


Program AUT CMD (SVF) 

CCR 

STH CPD 


Unsigned 12-bit integer results will be stored in the 1st, 2nd or 3rd word of each 3-word 
block as labeled above. 
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I:T:1 - l:T:2 - Output source is one of the 1024 word time interval data fields labeled 
'TRAIN 1 1 and 'TRAIN 2' respectively. If less than 1,024 spikes have been recorded 
(counted), non-required locations are set to zero and punched too'. (See program 
description IC1 and IC2'.) 

NOTE: - File designation terminating RETURN initiates the output subroutine to punch a 1 ft. 
leader tape. Therefore the Low Speed Punch (LSP) control should be turned ON previous 
to the RETURN instruction! 

NOTE: - Non-significant character printout during data punch output is due to parallel 
teletype/punch logic and is of no importance. 

At the end of data output, a 1 ft. trailer tape is punched. The LSP control should then 
be turned OFF, otherwise all subsequent teletype communication is punched on paper tape! 

Program control initiates a new "OPT:" printout and returns to PRSi 

Each 12-bit word requires two consecutive punch columns. The first half word is punched 
on tracks 1 - 6 of the first, the second half word on the same tracks of the second column! 

NOTE: - MAKE SURE THAT EACH PAPER TAPE IS LABELED WITH SUFFICIENT DATAS 
FOR ADDITIONAL PROCESSING! 

4.1.5 L 

Code terminating RETURN initiates program request typeout "T". The user is expected 
to design one of the time interval data fields, referred to as 'TRAIN 1' and 'TRAIN 2' 
respectively, typing a "1 " or "2" and terminating RETURN. 

The program then evaluates five of the most required statistical measures of central 
tendency. 

N - Number of intervals (spikes) of the train designed. 

T - Sum of time intervals, i.e., total record duration, 

MU - Interval mean value. 

S - Standard deviation, 

C - Estimate of the coefficient of variation (having the value unity for an exponential 
population)! 





Formulas of computation 
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C = S / MU 


Output format - Data output is in standard decimal floating point format. 


+0.mmmmmmmE+ee 


Example 


I_exponent (base 10) 

—__..mantissa (decimal) 

(time base 0.001 sec) 


N +0.1013000E+04 

T +0.2840200E+05 

MU +0.280375] E+02 

S +0.1816645E+02 

C +0. 11 77067E+02 


1013 spikes have been counted 
over a period of 28,402 sec. 

Interval mean is 28,03751 msec 
with standard deviation 18,16 msec 
Estimate of the coefficient of 
variation is 11,77 and indicates 
strong rejection of the Poisson hypothesis'. 
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4/2 IC1 Interval Counter 1 channel 
IC2 Interval Counter 2 channel 

These two programs are AX08 interface monitor subroutines and perform digitalization of 
interspike time intervals in series of 1024 events (spikes) per channel. 

The programs introduce with label and input linkage mnemonics. 


AX08 Schmitt-Trigger W501 inputs SI and S2 accept negative going, negative spike pulses. 
The switching thresholds can be varied over the range of 0 to -2.5 V by front panel potentio¬ 
meters (see program OPT, option ST). Schmitt-Trigger input S3 accepts the counter start pulse 
which may be set either by an optional pulse device or by a synchron pulse track on the tape 
recorder device. In either case, the pulse specifications mentioned above should be 
recognized 1 . 


Start pulse enables the RC-clock and initiates counting procedure. Each event (spike) 
recorded will save the clock counter value on consecutive train storage locations as a 12-bit 
unsigned integer, that represents the RC-clock pulses counted during the last interval. No 
zero interval will occur 1 . 


NOTE: - IF ANY TIME INTERVAL EXCEEDS 4095 (multiplied by the time base!) AN 
OVERFLOW ERROR MESSAGE WILL OCCUR'. 


To avoid this interrupt make sure that the RC-clock calibration allows counting of the longest 
interval recorded 1 (See 4.1.1) 

IC1 counting procedure terminates: 

- when 1,024 or N 1,024 (see 4.1.3) time intervals have been digitalized and stored. 


- when 1024 or N intervals have been processed by one of the two counters in IC2 mode. 

- after OVERFLOW message printout. 


- when a data break has occurred, i.e., when the tape recorder playback has been turned 
off. (This will cause an OVERFLOW error message because the time period following the 
break point will be interpreted as an interval which exceeds 4095'.) 

Non-required storage locations are set to the value zero and the trains then saved on a full 
(or half) DF32 disk track. 

Disk transfer completed will cause an end-of-job mnemonic typeout, and program control turns 
on PRS to receive the next program code. 

The user is recommended to proceed then with program OPT, option L, to get general 
information about raw data. This will be helpful to plan additional analysis! 


9 





4.3 Autocorrelogram 

This program computes the order-dependent statistical measure involving time intervals 
between non-successive events, known as the estimates of the renewal density function or 
autocorrelation function. 


Principle of operation: The range of observed interval lengths (see X SCALE FACTOR) is 
divided into 128 bins of equal class size d; if the i-th observed interspike interval X. satisfies 
the inequalities 


(j-l).d < X. < j.d 

then that interval is placed in bin j of the histogram. The bins are numbered 1, 2,.. 
For higher-order density estimates, the inequalities hold 
(j-l).d < X?< j.d 


128. 



where 



r 


i+n-l 

k=i 


X, 


1, N-n+1 


and the estimate is of order n! 

So n histograms have to be evaluated and summed up. To compute the average range of a n-order 
autocorrelation estimate, use the rule 


T O' n.MU 
n 


(for MU, see program OPT, option L) 


Since the renewal density is the sum of_a_M orders, 
decrease for increasing n*. 


i.e., n_j.o= , the estimate error will 


Contradictory to this rule is the fact that higher order histograms will flatten out due to the 
finite number N of the interval sample. 

So, for example, if you take n = N, the n-th histogram will be only one bin (the last one) 
incremented, which is of less weighting value for the overall estimate! 



n=l\|/2 


n=N 

-^Ss. 


(The figures show the k-th order histograms for k = n) 


Use - Code terminating RETURN initiates label and request printout. 
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TRAIN - The user is expected to design the spike train to be processed (Typing 1 or 2) and 
terminating RETURN 1 . 

LAG ORDER - The user is expected to type the order n of the autocorrelation function dis¬ 
cussed above. 

LAG ORDER 1 initiates computation of the first-order autocorrelogram estimate, known as 
"interval histogram 1 ". 


X SCALE FACTOR - The user has to specify the time of observation necessary for analysis, 
according to the table listed below: (TIME BASE 1 mSEC) 


X SCALE FACTOR 

CLASS WIDTH 

TIME OF OBSERVATION 

0 

1 mSEC 

128 mSEC 

1 

2 

256 

2 

4 

512 

3 

8 

1024 

4 

16 

2048 

5 

32 

4096 

6 

64 

8192 


Terminating RETURN then initiates processing of n histograms (see principle of operation) 
which are summed up in the first words of a 1 28 x 3 integer array (see REGISTER MEMORY 
CHART, 4.1 .4). After processing completed, AUT calls up the histogram display subroutine 
HST (see 4.4). 

The user may become troubled when he will determine the order n of an autocorrelation 
analysis. As dsmonstrated above, each order will shift the histogram to the right with a 
delay of about the mean interval time MU. 




If n = 5, the user will have to expect an estimate error (which may be small at all) due to the 
classes lost of 6, 7 or even 8 order histograms'. To avoid this error, start the analysis with the 
n calculated from the formula given in the section "principle of operation', then recall AUT 
and increase n until the autocorrelation estimate will show no significant variation. 
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Analysis - The autocorrelation is a probability density function and specifies the probability of 
encountering any event (spike) as a function of time after a given event (spike). 


The first-order autocorrelation function, i.e., the "interval histogram" specifies the probability 
of encountering the next event (spike) as a function of time after a given event (spike)'. 

If the series of spikes is the realization of a Poisson process, the autocorrelogram will flatten 
out to the value 

7 = mD - ” 

For detailed information see REFERENCES APPENDIX in this manual 1 . 

4.4 HST - Histogram Display Subroutine 

HST may be called in memory either by the user or autonomously by programs AUT, CCR, 
and STH.HST introduces with label - and 

Y SCALE FACTOR 

request printout. The user is expected to specify the vertical counting range for histogram 
class display. He may use the following table: 


Y SCALE FACTOR 

Y RANGE 

0 

0- 127 

1 

0- 255 

2 

0- 511 

3 

0- 1023 

4 

0- 2047 

5 

0- 4095 


The program then displays 128 classes of counts, the class width pre-selected by programs AUT 
CCR, and STH. 

If one or more classes reach the top of the display area, which is marked by half- and full-scale 
dots, the user should recall HST, which is in PRS state, and then increase the Y SCALE FACTOR 
by one. The display will then be scaled by a factor 2 

NOTE: - THIS PROCEDURE MAY BE REPEATED UNTIL NO CLASS EXCEEDS THE UPPER 
SCALE RANGE LIMIT'. 

4.5 CMD - Cumulative Distribution Function 

Principle of operation - The program computes and displays the estimate of the probability 
distribution 
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t 

F(t) = £ f(T).dT 
o 

where f(T) is the probability density function discussed in 4.3'. 
The program evaluates 


F. 

i 


Tn XI n. 

n y | 


i = 1,128 


from the class frequencies n. stored in the first column locations of the register array (4.1 .4) 
and deposits F. in the 1 second column locations of the register. CMD will therefore 

make referenc4 to programs already evaluated, as AUT (LAG ORDER 1) or STH. 


Use - Code terminating RETURN is immediately followed by CRT display and a vertical scale 
grid, marked from 0 to 1 by steps of 0.051 

NOTE: make sure that 'sth 1 or 'aut' 9lag order 1) HAVE BEEN PROCESSED BEFORE USE 
OF "CMD" ’. 


Analysis - The estimates of the probability distribution is known too as "cumulative frequency 
distribution", and measures the probability that an event (spike) will have occurred by time t 
from the last event (spike). Then, most frequently, the equality listed below is used: 


p(T 1 ^T<T 2 )= FOy-F^) 

It will give the probability that the interspikes interval time will be in limits T^ and 
respectively. 

4.6 CPD - Conditioned Probability Density 


Principles of operation - CPD computes and displays an estimate of the order-independent 
conditioned probability (postpulse probability, hazard function). 


'f (t) 


_ f(t) 


l-F(t) 


The program evaluates 


i 


N ■ 2_ n. 
sr 


‘AT 


i = 1,128 


when 


n. = interval class frequency 

N = number of intervals counted 
AT = interval class width 
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Use - Code terminating RETURN is immediately followed by CRT display and the grid, 
described in 4.5'. 


NOTE: - MAKE SURE THAT "AUT" (LAG ORDER 1) OR "STR" HAVE BEEN PROCESSED 
BEFORE USE OF "CFD"'. 


Analysis - The condition density function specifies the probability density of encountering the 
next event (spike) at the time t after the last event (spike), given that there was no event prior 
to the time tl 


For a Poisson process, the condition density function is a constant. 


* w - * ■ md 


4.7 SCG - Serial Correlogram 

Principle of operation - The serial correlogram is the set of the serial correlation coefficients 


R ; = £«x.-/*->(x. +j -yw-» (-..,01,0,1 


SCG computes 128 coefficient estimates 



N-i £ X i X i+j ~ (N-i) 2 [ £ X j) ( SI X i+j) 


( NFT^ X f" 


(TFj ) 2 (^ X i) ) 


- 7m^ 2 (XX.^. 

N-| i+| (N-| i+| 


for j = 1,..., 128 


In this subprogram, all computation is performed in three-word floating-point mode. The 
coefficients are held in a 128-floating-point-array (see REGISTER MEMORY CHART, 4.1. 


4) 



Use - Code terminating RETURN initiates label printout and request 
TRAIN 


The user is expected to design the spike train to be processed, typing a 1 or a 21 

The program is now running inOLS state for about 5MIN (job time depends on actual data 
file length Nl) Processing complete, control turns on URS and initiates CRT display. 128 
correlation coefficients are displayed, together with a scale grid, marked from -1 to +1 by 
steps of 0.11 

Analysis - Serial correlation coefficient analysis is explained with respect to the well-known 
joint interval density . 
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Joint interval times are displayed in the form of a scatter diagram, in which the length of a 
joint-pair of intervals is represented by abscissa and ordinate respectively. 


n 




x i 



Figure A) shows an interval histogram of a set of 27 intervals; the same distribution was used 
to reconstruct two sequences, with scatter diagrams E) and C) for the joint interval distribution 
sketched above. Intervals of train B arise from a renewal process and therefore are drawn 
independently from the distribution A), the scatter diagram will look quite "random. 11 

The intervals of train C arise from the same distribution, but have been arranged in a decreasing 
sequence; the related scatter diagram shows a "non-random" pattern. 

The problem arises from 'subjective' pattern recognition and interpretation of more or less 
significant patterns, i.e., correlations. A quantitative measure of such a correlation is fur¬ 
nished by the serial correlation coefficient defined above, and we will illustrate the case for 
R,, i.e., the joint interval serial correlation (for nonadjacent intervals the conclusions will be 

To illustrate the correlation product 

c = (X. - MU) (X. +1 - MU) 

dotted lines on scatter diagrams mark mean interval times 
MU on both abscissa and ordinate and divide the diagrams' 
first quarter in four sets of scatter points (i.e., joint interval 
pairs). 

The related correlation products will be positive or negative 
signed according to the set, the joint pair (X.,X. + ^) belongs to: 


similar). 



For the sequence C this means, that almost all (except one) products are positive, and the 
correlation coefficient will not be far from the standard deviations square value. 

For the independent sequence B there will be almost the same number of positive and negative 
products to be summed up, and the coefficient therefore to close to the value zero . 
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The expected value of the serial correlation coefficients of all orders (1 - 128) is approxi¬ 
mately zero if the intervals are drawn independently from a common distribution. 

Positive contributions to serial correlation coefficients are due to long-term trend in data 

For detailed analysis see REFERENCE APPENDIX. 

4.8 CCR - Cross-Correlogram (Cross Correlation Histogram) 

Principle of operation - The cross-correlation histogram is a measure to test whether two spike 
trains are independent. 


Train I is the reference train or, in terms of the 'poststimulus time histogram 1 , the "stimulus 
train. " Forward and backward recurrence times are then counted in train 11 relative to any 
event (spike), or ' stimulus, 1 in train I .A histogram will be computed for all spikes in train 1 
to all spikes in train II 



Use - Code terminating RETURN initiates label and request printout. 

X SCALE RACTOR - 

The user has to specify the time of observation. Forward and backward recurrence times 
will be counted, so the origin of the time scale will be between class NO. 64 and NO.65. 


(TIME BASE 1 mSEC) 


X SCALE FACTOR 

CLASS WIDTH 

TIME OF OBSERVATION 

0 

1 mSEC 

- 64 mSEC — 

+ 64 mSEC 

1 

2 

- 128 — 

+ 128 

2 

4 

- 256 — 

+ 256 

3 

8 

- 512 — 

+ 512 

4 

16 

- 1024 — 

+ 1024 

• • * 

# • • 

• •it 

• • • • 


Analysis - The cross-correlation histogram is an estimate of the cross-correlation function 
defined as the.probability density of encountering any event (spike) in train 2 as a 
function of time after a given event (spike) in train 1. 
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The cross-correlation histograms have two uses: 


- as a corroboration of independence, indicated by a flat cross-correlogram. 

- as a means of exploring suspected interactions, indicated by positive or negative going 
elevations (peaks), with respect to the reference level. 

The following symmetry relationship for the cross-correlation holds: 

X 1 2 (0 ("0 

MUj = MU^ 


where Ml^ and Ml^ are the mean interspike intervals of train 1 and 2 respectively. For this 

reason, the cross-correlation estimate is measured in only one direction, from train 1 to train 2, 
but for both positive and negative values of t. 

Prediction of correlation or independence should be carefully examined; errors arise from 
different sources, some of them being mentioned below. 

- independently firing pacemakers of near the same frequency exhibit peaks in the cross¬ 
correlation histogram, which oscillate with the period of the pacemaker neurons'. 


- distinction between interaction (direct or indirect) and common input. 



- weak dependence will be indistinguishable from sample 'noise". 

- combinations of positive and negative interactions may coincidentally combine and evoke 
flat cross-correlograms, i.e., prediction of independence will be guessedl 

For further and detailed analysis, see REFERENCE APPENDIX '. 

4.9 TIS - Time Interval Sequence 

Principle of operation - Most of the statistical analysis done will deny the 'null hypothesis', 
but rather leads to the prediction of 'dependence', 'functional relationship' or 'significant 
patterns'. Although these measures will do so, they will not give any information about the 
source of the relationship or the pattern itself. This information is lost during computation 
and up to the time now, no ' all- round' pattern-recognition program exists'. 

In practice, most of the prediction concerned with these questions are drawn out from raw data, 
i.e., from the original sequence of events (spikes). 
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To help the user, this program will display the interval times between adjacent events as a 
function of their number of occurrences in the sequence of events. This is done on a small 
area, compared with any other record representation, and the intervals are shown as ordinate 
values, a more accurate form than any numerical listing. 

Use - Code terminating RETURN will cause label and request printout: 

TRAIN 


The user is expected to specify the sequences of events (train) to be displayed. Control 
immediately turns on URS, and the sequence of intervals appears as a set of dots on the CRT 
device. 


The sequence and the vertical scale factor may be controlled on-line by the use of ADC 
control knobs on the AXJ08 front panel as follows: 


ADC 1 - The first interval, out of the sequence to be displayed, may be specified, its 
number of occurrence appears in the display record field "C" as an even integer in the 
limits of 0 to 1022. 


ADC 2 - Up to 512 intervals may be displayed simultaneously. This upper limit may be 
scaled by powers of 2, down to the smallest set of 4 consecutive intervals. The specified 
number appears in the field "X" of the display record. 


ADC 3 - Vertical timescaling by powers of 2 is performed by the use of the third control 
knob. Realtime values may be calculated with respect to the time base, selected during data 
conversion (see OPT, 4.1.1). Fullrange numbers from 64 up to 4096 time units are optional 
and appear in field "Y" of the display record. 


NOTE: - THE SUM OF THE TWO CONTROL COUNTERS "C " AND "X" SHOULD BE 
CHECKED TO HAVE A VALUE LESS THAN THE TOTAL NUMBER OF EVENTS (SPIKES) 
COUNTED FOR THE TRAIN SPECIFIED (N^1024). 


Otherwise, non-significant CPU core data will be interpreted and displayed as time 
intervals too' 


4.10 STH - Stimulus-related Time Histogram (Post-stimulus time histogram) 

Principle of operation - Forward recurrence times are counted with respect to each stimulus 
event. The recurrence times then are re-arranged in a histogram. This post-stimulus time 
histogram is an estimate of the cross-correlation function described in section 4.8, but with 
the restriction that calculations are performed for positive times and in one direction only. 

Use - Code terminating RETURN initiates label and request printout; 

X SCALE FACTOR 

The user is expected to specify the time of observation, according to the table listed in 
section 4.31 
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Request terminating RETURN initiates data linkage mnemonics printout: 

51 stimulus events data channel 

52 stimulus-related events data channel 

53 start pulse channel 

NOTE: - STH is an on-line analysis subroutine and will not process interval sequences 
counted by the use of program IC2 1 . 

The counter stops, when either 1024 events from input channel S2 or N stimulus events from 
input SI have been recorded. The number of stimuli (N) must be specified previous to the 
use of STH (see program OPT, 4.1.3). 

With increasing time of observation, more events are counted after a given stimulus; though, 
W by the limit of 1,024 (which will be most frequently reached before N !) the number of 
stimulus events counted will decrease. 



The program autonomously calls subroutine HST for histogram display, when the counting 
procedure has been terminated. 


Although the program was designed for on-line analysis, the post-stimulus times counted will 
be saved as a disk file and may be called in core memory for computation of the measures of 
central tendency (program OPT 4.1.5) later on. 


Analysis - As mentioned before, the STH histogram is an estimate of the probability density 
function equivalent to the cross-correlation function between the train of stimuli and the 
train of spikes. If the stimulus has no effect, the histogram will be flat. 


5. CRT Display and Photo Recording 

The package was designed for on-line application and therefore high priority was given to 
the CRT information display concept. 

Each analysis program outputs data in the form of a CRT graphic, always composed of: 

- a record (head) 


- a scale grid 
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and a function or histogram plot 


The scale grids and the functions have been mentioned and described in the section 
PROGRAM DESCRIPTION. 


The display record is composed of a set of mnemonics and numbers, and was designed for 
photo record identification. 


0 

TIS T:0002 

C:0098 


X:0256 
Y:0128 


© 


STAP-8 

30.5.69 



1 . Program mnemonic refers to the program who's output is displayed. 


2. Package mnemonic 


3. Comment record field: 10 characters may be displayed in this field for additional 
user comments. During CRT display mode, program control always is in URS state, and 
therefore the teletype is sensed for comment or program code characters. The first key stroked 
CRT display is interrupted until the terminating RETURN key has been sensed. 


4. Photo record counter: This two-digit decimal counter is incremented after the photo 
record sequence, described below, and may not be reset. 


T Spike train number (with reference to Schmitt-Trigger inputs SI and S2), specified as 
input data file for the running program. 


C The scratch constant has a different meaning, according to the program running: 


AUT 

Lag Order n 

CMD 

- 

CPD 

- 

SCG 

Number of coefficients computed 

CCR 

- 

TIS 

First interval (no. of occurrences) 

STH 

Number of stimuli counted 


X Indicates the full range of the abscissa and must be interpreted according to the running 
program either as a number or as a time (multiplied by the time unit) . 

Y Indicates the full range of the ordinate and must be interpreted either as the number of 
counts (histogram), a real value of a function (CMD, CPD, SCG), or a time as in program TISl 

NOTE: - IF THE INTEGER NUMBER OF EITHER THE RECORD FIELD "X" OR THE RECORD 
FIELD "Y" EXCEEDS THE VALUE 4095, A DATA OVERFLOW OCCURS AND THE CHARACTERS 
0000 ARE DISPLAYED'. 

In this case, the user should refer to the teletype record of the value missed. 
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5.1 Photo Record Sequence 


S3 - If during the display sequence a -3 V pulse is sensed at input channel S3, the program 
control interprets the pulse as a user request for remote photo control. This option was 
designed for two reasons. First, function and histogram display have different sweep times 
(CRT refreshment), and therefore need different shutter interval (exposure time). 

Second, there is a package version available, where all the histograms make use of a 
storage CRT device. This option is used to avoid CRT flicker, due to histogram refreshment 
intervals. 


The record sequence initiated, CRT remote control is set to NON-STORE MODE to obtain 
clear-cut photo records. At the end of the record sequence, control again turns on STORE 


MODE. 


XR 20 - The photo sequence exposure counter is incremented for each refreshment sweep. 
During this sequence, the AX08 output level option XR 0020 is set to the logic value 1 and 
cleared immediately after conclusion of the counter loop. This option may be used for photo 
remote control (shutter logic B). 


NOTE: - THE EXPOSURE COUNTER (STORAGE LOCATION 0322) MAY BE RESET 
ACCORDING TO THE CRT INTENSITY AND RECORD FILM SENSITIVITY AVAILABLE 1 . 
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Photo Sequence 
Flow Chart 



Function Display 


Histogram Display 
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APPENDIX B 


ERROR DIAGNOSTIC 


INCORRECT PROGRAM CODE 


? 


DF32 ERROR MESSAGE 


OVERFLOW AT 0000 


ZERO COUNT LOAD 


Error diagnostic messages may occur 
and are listed and discussed below. 

Loading procedure error messages are 
discussed in section ERROR DIAGNOSTIC 
of the PD P-8/1 console manual 1 


TT intercom error messages 

A character string with a leading was checked 
to be incompatible with valid program codes. 

State PRS 


A character string was checked to be composed 
of more than 16 characters allowed. 

State PRS 


Disk transfer error messages 


A DFSE skip request was missed during disk 
transfer. The error is due to hardware break¬ 
down, and the systems engineer should be contacted. 
State HLT (program flow interrupted) 


Program error messages 


The link bit was sensed for data overflow. 
The result of an arithmetic operation will 
not be valid. 

State program continuation 


IC1 at 3475 

IC2 at 3504 and 35]0 


(An interval counter has counted 4,096 time units) 
CMD at 3462 

(divisor smaller than dividend) 


A counter location was initialized and 
set to the value O. 

State PRS 

OPT at 3425 
AUT at 3405 
CCR at 3407 

A non-existent data file (spike train) 
had been specified for processing. 
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Program 


STAP-8 


Page 

1st A 

last A 

Content 

0 

0 

177 

CONSTANTS/SR POINTER/SEQ CODE LIST/SEQRDD 

1 

200 

377 

MESAGE/PHOTO/RDD/DFERR 

2 

400 

577 

MESSIN/OFLOW/OCTDUM 

3 

600 

777 

MASTER/BELl/BINDEC 

— 

4 

1000 

1177 

DECBIN/WRT 

5 

1200 

1377 

TRANSF/GRID1/BINDUM/LOADTR/WC/CA/SAVE/LAYOUT 

6 ^ 

1400 i 

1577 

LABEL/SAVETR 

7 

- ■ j 

1600 | 

1777 

ADISP/GRID2/ZERON 

10 

2000 ; 

2177 

FLPOUT/READC 

11 

2200 

2377 

CHDSC 

— 

12 

2400; 

2577 

DECTBL 

13 

2600 | 

2777 

(BINARY PUNCH SA:2665) 

14 

3000 

3177 


15 

3200 

3377 


, 1 

i 16 

1 . 1 1 M | 

3400 

3577 

NON-RESIDENT LOADER/DF32 SAV^SEQUENCE/SA:34#0 

' 17 

3600 

3777 

OPT 1 

20 

4000 

4177 

OPT II 

21 

4200 i 

4377 

IC1 

* -1 

22 

,, , , ... .. . L— 

4400 1 

4577 

IC2 

■ 23 

4600 

4777 

HST 

24 | 

5000 

j 

5177 

AUT 1 

25 

5200 

5377 

j AUT II 

26 

5400 

5577 

CMD 

- 

; 27 

5600 

5777 

| TIS 

_30 

6000 

6177 

CPD 

! 31 

6200 

6377 

SCG 1 

32 

6400 

6577 

SCG II 

; 33 

6600 

6777 

I CCR 1 

34 

7000 

7177 

CCR II 

35 

7200 

7377 

STH 1 

36 

7400 

7577 

STH II 

37 

7600 

7777 

SAVE AREA FOR RIM/BIN LOADER/MONITOR RESIDENT 
i 


Remarks 


LOAD STEP I 
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Program 


STAP-8 


Page 

1st A 

last A 

Content 

0 

0 

177 

ZERO-PAGE 

1 

200 

377 

SR 

2 

400 

577 

SR 

3 

600 

777 

SR 

4 

1000 

1177 

SR 

5 

1200 

1377 

SR 

6 

1400 

1577 

SR 

7 

1600 

1777 

SR 

! io 

2000 

2177 

SR 

LLL , 

2200 

2377 

CHDSC 

i 12 

2400 

2577 

DSCTBL 

13 

2600 

2777 

FLPBUF (128 3-WORD FLOATING POINT REGISTER) 

14 

3000 

3177 

FLPBUF 

Lll. 

3200 

3377 

FLPBUF 

16 

3400 

3577 

RESIDENT SEQUENCE STORAGE LOCATION 

__ 17 

3600 

3777 


20 

4000 

4177 


Lll 

4200 

4377 


22 

4400 

4577 


23 

4600 

4777 


24 

5000 

5177 


25 

5200 

5377 


26 

5400 

5577 


27 

5600 

5777 

8-5-S/8-25-F FLOATING POINT PACKAGE 

30 

6000 

6177 

8-5-S 

31 

6200 

6377 

8-5-S 

32 

6400 

6577 

8-5-S 

33 

6600 

6777 

8-5-S 

34 

7000 

7177 

8-5-S 

35 

7200 

7377 

8-5-S 

36 

7400 

7577 

8-5-S 

37 

7600 

7777 

SAVE AREA FOR RIM/BIN LOADER 


Remarks LOAD STEP tl 
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Program 


STAP-8 


Page 

1st A 

last A 

Content 

0 

0 

177 

ZEROPAGE 

1 

200 

377 

SR 

2 

400 

577 

SR 

3 

600 

777 

SR 

4 

1000 

1177 

SR 

5 

1200 

1377 

SR 

6 

1400 

1577 

SR 

7 

1600 

1777 

SR 

10 

2000 

2177 

SR 

11 

2200 

2377 

CHDSC 

12 

2400 

2577 

DSCTBL 

13 

2600 

2777 

FLPBUF 

14 

3000 

3177 

FLPBUF 

15 

3200 

3377 

FLPBUF 

16 

3400 

3577 

RESIDENT SEQUENCE STORAGE LOCATION 

17 

3600 

3777 

TRAIN 1 

20 

4000 

4177 

/ .. 

21 

4200 

4377 

/ 

22 

4400 

4577 

/ 

► 23 

4600 

4777 

/ 

24 

5000 

5177 

/ 

25 

5200 

5377 

/ 

26 

5400 

5577 

/ 

27 

5600 

5777 

TRAIN 2 

30 

6000 

6177 

/ 

31 

6200 

6377 

/ 

32 

6400 

6577 

/ 

33 

6600 

6777 

/ 

34 

7000 

7177 

/ 

35 

7200 

7377 

/ 

36 

7400 

7577 

/ 

37 

7600 

7777 

SAVE AREA FOR RIM/BIN LOADER 


INTERVAL COUNTER MODE IC1 IC2 CCR 
27 


Remarks 






















































Program STAP-8 


r~ 

Page 

1st A 

last A 

Content 

0 

0 

177 

ZEROPAGE 

> 

200 

377 

SR 

2 

400 

577 

SR 

3 

600 

777 

SR 

1- 

4 

1000 

1177 

SR 

5 

1200 

1377 

SR 

6 

1400 

1577 

SR 

7 

1600 

1777 

SR 

10 

2000 

2177 

SR 

11 

2200 

2377 

CHDSC 

12 

2400 

2577 

DSCTBL 

13 

2600 

2777 

FLPBUF 

14 

3000 

3177 

FLPBUF 

15 

3200 

3377 

FLPBUF 

16 

3400 

3577 

RESIDENT SEQUENCE STORAGE LOCATION 

17 

3600 

3777 

TRAIN 1/2 

20 

4000 

4177 

/ 

21 

4200 

4377 

/ 

22 

4400 

4577 

/ 

23 

4600 

4777 

: / 

24 

5000 

5177 

1 / 

25 

5200 

5377 

/ 

26 

5400 

5577 

/ 

27 

5600 

5777 

8-5-S/ 8-25-F 

30 

6000 

6177 

/ 

31 

6200 

6377 

: / .. . 

32 

6400 

6577 

/ .... 

| 33 

6600 

6777 

/ 

34 

7000 

7177 

/ .. 

35 

7200 

7377 

/ ....... 

36 

7400 

7577 

/ 

37 

7600 

7777 

SAVE AREA FOR RIM/BIN LOADER 


Remarks SINGLE TRAIN ANALYSIS / FLOATING POINT MODE 


V 
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